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A  Big  Subsonic  Channel  as  a  part  of  the  U-15T-2  set  is  intended  to  investigate  the 
behavior  of  heat  protection  large-scale  samples  (~300  mm  in  length)  in  the  environment  of 
high-enthalpy  gas  flow  at  turbulent  boundary  layer  conditions. 

The  test  scheme  is  shown  in  Fig.  1.  The  sample  has  a  shape  of  spherically  blunted  long 
cone  with  a  cover  made  from  a  studied  material.  The  sample  is  placed  into  a  cone  shroud 
connected  with  the  heater.  The  sample  is  fixed  on  a  rod  with  a  strain-gauge  balance,  the  rod  is 
connected  with  a  device  of  the  sample  movement. 


Fig.l.  BSC  test  scheme. 

1  -heater;  2  -  channel;  3-  tested  sample;  4  -  strain-gauge  balance;  5  -  back  protection  shell; 

6  -  sample  moving  device. 

The  channel  design  provides  that  the  area  of  the  gap  between  the  model  and  the  channel 
cones  is  closely  the  same  along  the  surface  and  sonic  flow  is  located  at  the  model  bottom.  The 
model  bottom  diameter  decreases  due  to  ablation  and,  to  keep  the  gap  area,  the  model  is 
automatically  moved  upstream,  the  movement  is  controlled  by  the  heater  pressure. 

When  heat  protection  cover  is  heated,  recession  rates  of  its  binder  and  reinforcement 
(fabric  as  a  rule)  are  different,  so  a  relief  develops  on  the  surface,  which  can  induce  additional 
forces  and  moments  due  to  interaction  with  the  flow.  The  aerodynamic  six-component  strain- 
gauge  balance  mounted  inside  the  model  is  intended  for  measurement  of  these  forces  and 
moments.  The  balance  design  particular  features  are  considered  in  report  [1],  the  main 
performance  characteristics  are  given  in  the  table  below. 


Characteristics 

Components 

X,  kN 

7,  kN 

Z,  kN 

Mx,  kN-cm 

My,  kN-cm 

Mz,  kN-cm 

Nominal  load 

30 

5 

5 

1 

35 

35 

Resolution 

0.05 

0.008 

0.008 

0,0008 

0.05 

0.05 

Main  error ,  % 

0,2% 

Typical  plots  of  forces  and  moments  measured  by  the  balance  are  presented  in  Fig.  2. 
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Fig.  2.  Variation  of  forces  and  moments  during  a  test. 

Surface  condition  of  tested  sample  is  studied  by  means  of  roughness  measurement  both  in 
longitudinal  and  transversal  directions,  measurement  of  shape  with  the  use  of  coordinate¬ 
measuring  machine  UPCM  (Opton  Co,  Carl  Zeiss,  Germany)  and  photo  filming  of  the  sample 
and  its  fragments. 

Influence  of  surface  roughness  on,  for  example,  rolling  moment  factor 
mx  =Mx/SL(pV2/2)  is  shown  in  Fig.  3.  Values  (p V2/l)  correspond  to  the  model  bottom,  S  and  L 
are  the  bottom  area  and  length  of  the  sample;  roughness  parameter  Ra  (average  profile  height  [2]) 
is  a  medial  one  between  that  in  longitudinal  and  in  transversal  directions.  Though  there  is  no 
obvious  correlation  between  factor  of  moment  mx  and  parameter  Ra,  some  tendency  can  be  seen. 
Obviously,  it  is  not  roughness  elements  height  itself  that  defines  aerodynamic  loads  but  orientation 
of  the  elements  coupling  with  orientation  of  the  model  in  the  channel.  If  roughness  elements 
orientation  is  random,  it  is  only  an  average  level  of  the  induced  loads  that  can  be  assessed. 

The  model  shape  measurement  gives  a  possibility  to  find  a  variation  along  the  sample  gen¬ 
eratrix  of  material  loss  thickness  5  calculated  as  a  distance  between  initial  and  final  shapes.  Figure 
4  demonstrates  a  typical  curve  5(s),  where  s  -  distance  from  the  stagnation  point  along  the  surface. 
Fig.4  shows  that  heat  flux  distribution  along  the  model  generatrix  is  essentially  non-uniform.  In 
particular,  the  sharp  increase  in  8  at  ,y~55-M50  mm  indicates  boundary  layer  transition. 

Based  on  the  after-test  shapes  analysis  and  corresponding  transition  criteria,  algorithm  for 
calculation  of  heat  transfer  in  the  subsonic  channel  had  been  developed.  It  was  assumed  that 
boundary  layers  of  the  sample  and  the  channel  develop  independently  till  the  point  where  the 
channel  boundary  layer  achieves  the  sample  surface  or  the  sample  boundary  layer  connects  with 
the  channel  wall.  Downstream  this  point  the  flow  is  considered  as  the  gap  flow  and  heat  transfer 
is  defined  by  mass  overage  parameters.  The  following  correlations  using  the  effective  length 
method  [3]  were  involved  for  laminar  and  turbulent  heat  transfer: 


Fig.  3.  Rolling  moment  factor  dependence  on  roughness. 
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Fig.  4.  Variation  of  material  ablation  thickness  and  heat  transfer  factor  along  a  sample  surface 
qwMm  =  0.332Pr  "2 /3  •  p  *  U  5  ■  Re  ,  *-°'5  {h  eJam  -Hw)-K, 


qw,turb  =  0.013  Pr-3/4p*[/5 -Re  a”  4  ~  » w) 

,  p*t/5/  PgC/gA  )RlP*u&*(He  -njfdx 

where  Re  */  = - — ,  Re  *  = - ,  /  =  ° - , 

P*  Pa  R2p*Up*{He-Hwf 

Kt  -coefficient  accounting  for  pressure  gradient  and  temperature  factor  ,  A  -energy  loss 
thickness. 

To  find  pressure  distribution  along  a  gap  between  sample  and  channel,  area  of  the  gap 
cross-section  is  calculated  as  a  tore  segment  and  pressure  is  calculated  by  iterations  to  fit  a 
given  mass  flow  rate.  Pressure  along  the  sphere  bluntness  is  calculated  as  P(x)=P0o(l-ax2+bx4), 
a  and  b  factors  are  determined  from  the  term  of  the  function  and  its  gradient  continuity. 

In  response  to  the  U-15T-2  jet  non-uniformity,  method  of  mass-average  values  ([4])  was 
involved.  Enthalpy  profile  H(R)  is  assumed  as 


H~Hn 

IT  _  IT 

11  axis  11  n 


=  1-1- 


-(l-R)iV/2  2, 


where  values  Hn  and  N  are  determined  from  the  term  of  mass  and  energy  conservation  lows.  In 
accordance  with  the  mass-average  values  method,  a  variable  of  flow  function  is  introduced  as 


1  R 

—  \Rp -U-dR 


and  dependence  of  total  enthalpy  on  flow  function  is  calculated  both  for  the  sample  and  the 
channel  as 


1  v 

H0=-\H{y)dy. 

Ilf  J 
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Fig.  5.  Total  enthalpy  dependence  Fig.  6.  Flow  parameters  variation 

on  flow  function  along  the  sample  surface 

Typical  dependence  of  total  enthalpy  on  flow  function  for  the  sample  and  the  channel  is 
shown  in  Fig.  5.  Calculations  were  conducted  for  the  case  of  pressure  in  the  heater 
/>oo=l-25MPa,  mass-average  total  enthalpy  H0 0=4.6  MJ/kg,  mass  flow  rate  G=1.8  kg/s.  As  an 
illustration  of  the  calculation  results,  Fig.  6  shows  variation  of  pressure  P,  flow  velocity  V  and 
density  p  along  the  sample  generatrix. 

The  calculated  curve  of  heat  transfer  factor  (a/cp)o  for  smooth  impenetrable  surface  is 
shown  at  the  plot  Fig.4.  Length  of  laminar  zone  is  defined  by  a  well-known  correlation 

Reg  =  A-e°'2M,  where  factor  A  is  chosen  equal  to  165  to  fit  after-test  shape.  Extension  of 
transition  zone  follows  from  correlation  Re  §  en/j  =  1 ,5Rc  §  [5],  It  can  be  seen  that  the 

assumed  mathematical  model  reflects  in  general  the  real  heat  transfer  variation  along  the  model 
generatrix. 

Experimental  data  on  measurement  of  heat  flux  and  pressure  variations  along  the  sample 
cone  surface  are  given  in  Figs.  7  and  8.  Figure  7  shows  variation  with  time  of  heat  flux 


Fig.  7.  Flcat  flux  hystory  measured  and  Fig.  8.  Pressure  variation  along  the  cone  surface, 

calculated  in  4  points  on  cone  surface  calculatted  and  measured 
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measured  in  4  points  along  the  surface.  Calculation  results  for  the  test  conditions  are  presented 
as  a  band  between  two  dotted  lines  on  the  same  plot.  The  upper  line  corresponds  to  the  point  5 
(50  mm  form  the  sample  bottom),  the  lower  one  corresponds  to  the  point  3  (the  cone  generatrix 
middle).  Correlation  of  experimental  and  calculated  data  on  pressure  distribution  is  given  in 
Fig.  8.  It  can  be  seen  that  there  is  a  rather  good  coincidence  of  experimental  and  calculated  data. 

The  presented  experimental  technique  for  Big  Subsonic  Channel  supported  by  the  flow 
parameters  and  heat  transfer  calculation  software  increases  essentially  possibilities  of 
experimental  investigation  of  heat  protection  material  behavior  in  turbulent  heat  transfer 
condition  and  of  ablation  surface  relief  influence  on  force  loads. 
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